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The effect of external colour field impact on the instanton liquid is studied. In the course of this 
study the corresponding effective Lagrangians are derived for both regimes of weak and strong 
external field and in long wave-length approximation. The example of Euclidean colour point- 
I like source is analyzed in detail and the feedback of field on the instanton liquid is estimated as 

' a function of source intensity. 

(N 
o 

Q I The declarations of discovering new state(s) of matter in relativistic heavy ion collisions at RHIC 
^ ■ which are actively wandering in the papers nowadays are sometimes based on the results of differ- 
ent nature. From one side it is the striking result of direct experimental measurements of a strong 
suppression (comparing to pp and pA collisions) of particle production at high transverse momen- 
, tum well-known as a jet quenching. And although the jet reconstruction in these experiments is a 
^ ^ nontrivial task the (and accompanying) result (s) is(are) interpreted as a degradation of hard parton 
Q ■ (initiating a jet) energy induced by medium (new thermalized matter) produced in collision long be- 
^ fore hadronizing in the QCD vacuum. On the other hand the convincing success of phenomenological 
^ : analysis of the other measurable characteristics based on the perfect liquid hydrodynamics results in 
^ the questions about the sort of plasma (if produced) and the origin of the QCD vacuum Ref. p. 
Ph| In the present paper we are trying to understand if the instanton liquid model of QCD vacuum 
o ' could be indicative in interpreting these RHIC results. The instanton liquid here is considered in 
^ ' the framework of the stochastic ensemble of instantons in the singular gauge and the generating 
functional is evaluated by the variational principle developed in the paper Ref. 0. Comparative 
• ^ . simplicity of the superposition ansatz and variational procedure allows us to analyze the effects 
^ I almost analytically. 

d • As a major configuration saturating the generating functional 



Z = Jd[A] e-^(^) 



where S{A) is a standard Yang-Mills action we take the approximate solution for the Yang-Mills 
equations in the form of the following superposition 



N 



A;ix)=B;ix)+Y.Alix;^,) , (1) 

i=l 

here implies the field of (anti-)instantons in the singular gauge 

Ki^) = - ^"^Vbt^u au{y) , av{y) = -T-, — 2 % ^ y = ^ - z ^ (2) 
^9 y +p y 

with the parameters ■ji = {pi, Zi, Ui) describing the i-th instanton of the p size centered at the pseudo- 
particle coordinate z, with the matrix of colour orientation u, and g denotes the coupling constant 
of non-abelian field. For the anti-instanton the 't Hooft symbols should be changed according to 



1 



f] ^ T], and the external field B'^{x) is viewed as well as the pseudo-particle field in quasiclassical 
approximation. 

The external field interaction with an individual pseudo-particle is defined by the following 
strength tensor 

= G%{B) + GIM) + GIM. B) , (3) 
where the first two terms correspond to standard strength tensors of non-abelian field 

GIM) = d,Al - d^Al + g r'^A^At , (4) 

with entirely anti-symmetric tensor f"''"^. In particular, 

GU^) = u;''%^^ M,^M,p WT7f ' 

where M^j, = J^j, — 2 y^y^, = j^^. The 'mixed' component of the instanton strength field looks 
like 

B) = g r'^BlAl - BlA^) = g r'^^uj^" ^ {B^ r^,,„ - B^ r^,,^) a„(y) . (6) 

Calculating now we receive the partial contributions of external field and each separate pseudo- 
particle as 

Gl. G% = G%{B) G;^{B) + G%{A) G%{A) + G%{A, B) G^M. B) + 

+ 2 GlAB) G%{A) + 2 G^AB) G%iA, B) + 2 G^M) G%iA, B) , ^ 

In order to keep the further steps as simple and transparent as possible we limit ourselves with 
the standard sum of partial contributions in the superposition ansatz action and hold the highest in 
IL density (precisely in packing fraction parameter np^) one particle contributions 

G?,., Gt, _ r . 



S{B,^)= jdx^^^^Y.jd^ 



The crossing terms of different pseudo-particles (which are proportional to the IL density squared) are 
neglected here because of very small packing fraction parameter characteristic to IL, i.e. np'^ ~ 0.01. 
Thus, the regularized generating functional for the IL model takes the following form (for denotations 
see Ref.|2]) 

r 1 /■ ^ 

^ = /^f^] M / n ^"'^''''^ • (9) 

The analysis can be easily done for the weak external field assuming the characteristic IL parameters, 
averaged pseudo-particle size p and the IL density n, unchanged and coinciding with their vacuum 
magnitudes. Those are fixed by some repulsive mechanism (see, however, the remark at the end 
of paper) for the particular choice of saturating configuration done above. In this case, from now 
on the integration over pseudo-particle sizes becomes unessential and all the pseudo-particles might 
be considered as having the same size p. It allows us to average over the pseudo-partile positions 
and colour orientations only by calculating the generating functional Q- Making use the cluster 
decomposition we obtain the corresponding average of exponential as 

(exp(-^))., = ( ? ^ ^^S'))^) ' (10) 
where {S,) = {{S^)), {S^S2) = {S^){S2) + {{S^S^)), . . .. 

The first cumulant is simply defined by the action averaged. Taking into account the direct form of 
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field strength tensors © and (jHl) it is evident that the following terms will only be present in the 
partial contribution after averaging over colour orientation 

(G^. = G%{B) G%{B)+{G%{A) G%{A))^+{G;M, B) G%{A, B))^+2 {G^M) G%{A, B))^ . 

Performing the colour averaging we use the equality 

{^'''^''') = , (11) 

implying as the number of colours. As a result we have that the last term 

16 

disappears due to the antisymmetric tensor f"-^'^, 

{G%{A) G%{A,B))^ = 0. 
Analyzing now the contribution of 'mixed' (repulsive) component 
G%{A,B) G%{A,B) = 4 r^^u'^" (Sj r^,,„ - B^ r^.^J a„(y) r'V^ {B^ r/„,, - B^, r^^,,) a,{y) , 
and averaging Eq. ljllll over the colour orientation (making use another equality) 

jabc jakc ^bk 

we have 

{GIM. B) G;:,(A, B))^ = (Bl - Bt nn,a) a^{y) (Sj r^,,, - B^ r^,,,) a,{y) . (12) 

Averaging over the pseudo-particle positions results in the following integral 

where 

/ = YP^ (13) 

because the basic IL parameters, as we agreed, are unchanged. Handling the 'mixed' component 
average we have it in the form as reads 



{GIM. B) GIM. B)U = y ^fl ^ Bl , (14) 



Finally, collecting all appropriate terms we find the effective action for the external field in IL as 

^^S)U.I,.(m^^^J^B')^Np, (15) 

here is the full number of particles in volume V with n = N/V and a single pseudo-particle action 
f] = Sn'^/g'^. The last term of Eq. ljlSp introduces the contribution of purely instanton component 
{G{A)G{A))^z- The contribution of repulsive term which fixes the pseudo-particle size in IL is omitted 
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in Eq. ()15j) so long as it is not a principal point in this context and adding it, leads to the insignificant 
correction to the last condensate term in Eq. . An amusing point is that the mass term of Eq. 
has been well-known for rather long time and as a matter of fact fixing the pseudo-particle size in 
the variational procedure of Ref.j2] is provided just by this mechanism of mass generation. With the 
characteristic IL parameters (iVc = 3 and number of fiavours Nf = 2) n/Ag^^ = 1.2, pAgcn = 0.27, 
(3 = 18 the mass estimate m ~ 0.44GeV looks pretty encouraging for p ~ IGeV and Aqcd in the 
interval of 200 — 300 GeV. The compatibility conditions for the equations resulting from Eq. ()15|1 is 
d^Bfj^ = which is satisfied by the pseudo-particle field Eq.© as well. 

Turning now to the next term of cluster decomposition to calculate the effective Lagrangian 
corrections we conclude immediately that in the second cumulant 

1 / / /" , G G /" , G2 G2 



2\\J d^^^jdx,- 



there are two nontrivial terms 



dx2 2 ""^ " " 1 , (18) 

here the index 2 underlines the fact that corresponding functions are dependent on X2- The remaining 
terms originate from either the interference terms (and are cancelled by the contribution of the first 
cumulant squared) or lead to the contributions anharmonic in B which are not in our interest 
for this paper. It was analyzed for the first time in Ref. [21 that G{B)G{A) in ((7j) generates the 
dipole interaction. However, this interaction does not manifest itself in the first term of cluster 
decomposition if the averaging over the colour orientation is performed. It comes into focus starting 
on the second order of decomposing. In particular Eq. ()17|l can be presented in the following form 



4 / 

2 J dx, dx2 ^^^^^^f"^^^^^ G'^M) , (19) 

if one exploits Eq.(jni) and Eq. (fTT|l keeping in mind that G|^j,(y4) 6*^^(742) is colour independent 
because of the identity uj°'^uj°''^ = 6^'^. Eq. lfTTH) should be also averaged over the pseudo-particle 
positions which results in the correlation function for the instantons in singular gauge developing the 
following form obtained in Ref. |lj 

where A = |xi — X2I, and for the anti-instanton the substitution e — > — e should be done. The 
analytical form of function Ig is not our priority here, however, it is shown in Fig.l. If the numbers 
of instantons and anti-instantons are balanced then the term proportional to the tensor e disappears 
and the contribution of correlator in IL is given by 

^ j dz G%{A) Glp{A2) = ^n (5^,5,^ - 5^p5,^) h . 

Now collecting the terms together we find the contribution of Eq. p7|) in the IL approach as 
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7n[ 



J dxidx2 Is (^^^ Gl.iB) Glp{B2) . 
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Figure 1: Correlation function Ig is given by solid line and the correlation functions Ji and J2 are 
given by the dashed lines. 

Clearly, it leads to an abatement of initial action and it is more convenient for analyzing to present 
the non-local factor of dielectrical susceptibility type in the Fourier components Ref. j3j 



dk 1 



16 1 



n Ukp) GUB{k)] GUBi- 



Numerical estimate of Is{kp) at the zero value of argument is 

and at Nc = 3, Nf = 2 the correction coefficient can be estimated as k 



^ n ls{0) 0.013. 

y c 



Analyzing now the term Eq. p8|) we present it as 



j dx2 2 



^0/3(^2) Gl^0(A2, B2) 



= ^ Q dx, dx2 GiM) viu, - VI,,) 



Qi-y X 



(20) 



and imply the dependence of G on the colour matrix u might be given by the common factor (without 
introducing new symbol for G). Formally, this term looks like the next one expanding in 1/Nc, i.e. 
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(~ cu^). However, using the identity for colour matrices 



fman , ,ak, ,nl ^kln , ,mg 
J UJ UJ = E Lk) 



we have 

a man , ,dk , ,nl rdbe , ,6c, ,ef\ 
UJ UJ f UJ UJ ') 

and then Eq. ()20|) receives the following form 



(21) 



j dx2 2 



dxidx2 \ Gl^{A) Glp{A2) Viu^Vips - G^^iA) G^fsW VkfSsViuj] a^a25 B'^B. 



m DTt 
2a ) 



The lower line here develops this form because of an asymmetric property of tensor G. Averaging 
over the pseudo-particle positions we may extract the correlation function in the following form 



dz 

V 

16 1 

7v 



G';,{A) G%{A2) niu^^nws - GlM) ^^^3(^2) nkmu^^ 



a^a25 



'A' 



^fia + J2 



'A' 



where A 



\X2—Xl\ 



is the unity vector. 



The simple algebra allows us to calculate the functions 



Ji= J dy 
dy 



J'. 



p' 


P' 


1 


1 


iy' + P'Y 




\y\ 


\z\ 


p' 


p' 


1 


1 


iy' + P'Y 




\y\ 


\z\ 



(16 - 8 t + 4 j9g + 6 (p^ + g^) t ~ 12 f pq) 
(4 + 5 t - 4 - 6 (j9^ + g^) t + 12 pq) , 



with z = y + A,t 



{z A) 



Similarly to Is we do not need their explicit forms here 



JyllJ'P M|A|'^ |^||A| 

but one may estimate their behaviours looking at the dashed lines in Fig.l. Finally, the additional 
contribution to the mass term reads as 



,,2^ 1 ~? ^ I dxidx2 



and in the Fourier components as 

^2 32 1 



'A" 



Sfia + J2 



'A' 
.P. 



A^A„ 



dk 



Bl{k)Bl{-k) 



(22) 



Estimating numerically the nonlocal correction to the mass we find out that Ji(0) ~ —1.4 un- 
likely above result. Then the mass term and corresponding correction in Eq. ()22j) come about at zero 
momentum Ott^ Nc and (— ^ 1-4), respectively. At the characteristic value /? ~ 18 it means the 
quantitative correction smallness or, globally, the corrections initiated by the second term of cumu- 
lant expansion are negligible at the contemporary values of basic IL parameters. There is another 
contribution to the effective Lagrangian which comes from the interaction of sources generating the 
external field with (anti-)instanton superposition 



TV 

Sint = E / c^a; j;(x)A;^(x;7i) . 
1=1 •' 
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Making use the cluster decomposition one expects the possibihty to calculate corresponding small 
contributions (if the sources are treated in the quasiclassical approximation) which are given by the 
correlation functions of the form {A'Jj^{x;'j)Al{y;j))^ Ref.|31. 

As the conclusion of this effective Lagrangian analysis Eq.()15p it is practical to address another 
approach to the problem Ref. [Sj. Let us suppose the quasi-classical field B is described in the infra-red 
momentum region by the initial Yang-Mills action without the term breaking down gauge symmetry 
as before. In particular, we consider the field of point-like Euclidean source of intensity e with only 
one non-zero n-th component 

5;(x) = (o,5"" v^ = |^^. 

Then B^ integrated over the 4-dimensional space gives 

/ dx (-^] = ^ X4 L , 
J \A7r \x\J An 

where X4, L are some formal upper limits of corresponding integrals. In this approach the contribu- 
tion of the first cumulant Eq. (fT3j) could be written down as 

{{S))^, = EX,, E = ^- + a L + l3nL\ (23) 

A n Tq 

with a = ^j^^zj 6^ np"^. The first term in defining E comes from the Coulomb energy of point-like 
source and tq represents a formal particle radius. The last term is originated by the gluon condensate 
and the previous term looks like negligibly small correction to the condensate term. However, this 
contribution linearly increasing with L is proportional to and has different physical meaning as 
a term additional to the self-energy of source. In other words, it demonstrates an impossibility for 
the source with an open colour to be available in IL because the amplitude of such a state is very 
strongly suppressed (e~'^) comparing to the condensate contribution if the screening effects are not 
taken into account. For the dipole in 'isosinglet' (s) and 'isotriplet' (t) states (i.e. Nc = 2) we obtain 

5»(x) = (0, s^'^),^ = ^ ir^i T ' 

^ An \\x — Zi\ \x — Z2\ J 

where Zi, Z2 are the dipole coordinates what leads to 

2 2 

( dx Bl = ^ X^l , I dx B^ = — X^iAL-l) , 
J An J An 

with I = \z\ — Z'2\ to be the distance separating sources. We have another confirmation of suppression 
effect for the states with open colour in IL, i.e. the energy of 'isosinglet' dipole state increases with 
/ enlarging and the corresponding coefficient is a ~ 0.6 GeV/fm if we take e ~ (7. 

Thus, we are quite allowed to conclude the regime of weak external field in IL is described by 
effective Lagrangian Eq. ()15p and basic IL parameters are within a well adapted interval. Moreover 
all the corrections originated by the second cumulant should be certainly neglected. 

Obviously, this conclusion will be considerably strengthened if a criterion of external field weak- 
ness is well defined and the point of how crucial is an assumption of the IL parameters unchanged is 
clarified. We are going to modify slightly the variational procedure of Ref . j2] to implement possibility 
of the changing IL parameters. We retain here the same designations to demonstrate precisely where 
the changes are introduced and imply S{B,'~f) in Eq.Q in the following form 

S{B, 7) = - ^ In d{p,) + 13 U,nt + E UeAli, B) + S{B) , (24) 
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The first term here describes one-instanton contributions with the following distribution function 
over the (anti-)instanton sizes 

d{p) = Cr^A^QCD (25) 
where b = ^N, - ^Nf, P = -Mn(AQCD p), 

4.66 exp(-1.68iVc) 



7r2(iV,-l)!(iV,-2)! 



The second term of Eq.(j23)) is responsible for providing pseudo-particles with repulsive interaction 
which fixes their sizes. The characteristic single instanton action is defined on the scale of average 
pseudo-particle size /3 = /5(p) where /3(p) = — IuCat^ — bln^AqcD p)- The partial pseudo-particle 
contributions grouped in the third term and we take only 



Uextili,B) = J dx 



because the other contributions at the standard IL parameters are small as we have seen. At last, 
the fourth term represents simply the Yang-Mills action of the B field 



S{B) = Jdx 
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The well-known property of exponential makes it possible to estimate the generating functional 
of Eq. Q with the approximating functional as 

Y>Y^ exp{-{S-S^)) , (26) 

where 

Yi = J D[B] TT^JUdl. e-^^(^'^)-^(^) , 5i(i?,7) = -Elnp(pO , 

i=l 

and yu(p) is an effective one-particle distribution function which may be derived with the variational 
procedure. In our particular situation a mean value of corresponding difference is given by 

1 1 „ 

N f 1 , ^ ^ M ^ f 1 1 rr r ^ r ^ r ^ 

^ J ^ d{/)^ 2 1? J f/mH7l,72) /i(Pl)/i(P2) + 



r N 1 r 

+ J '''' V jr, J '''' ° 

= j,.„ 1 + f „(^)%C7S^) , (27) 



with ( = wzji = T N^-i '^'^^ 1^0 = I dp p(p). Here we estimate the functional in the adiabatic 
(long wave-length) approximation. It means we consider the IL elements of some characteristic size 
(of the same order of magnitude as the mean distance between pseudo-particles) being equilibrated 
by the presence of some fixed field B. Then calculating the optimal configurations of pseudo-particles 
we found out the effective action in the mean field. Eq. ()27j] is given just in the form underlining 
that an integration is performed over liquid elements and the proper parameters describing their 
states could be dependent on the external field, i.e. could be the functions of coordinate x. Physical 
meaning of such a functional is quite transparent, it implies that each separate element of IL possesses 
a characteristic aptitude of screening external field assessed by Uext- 
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Figure 2: Pseudo-particle mean size as a function of applied external field. 

Calculating the variation of {S — Si) in /i(p) we have 

^(p) = C d{p) e-(<V+^ , 

where C is an arbitrary constant and we fix it demanding the coincidence of its value when the 
external field is absent with its vacuum average. Then 

M = C^J'^'^A'^cdp'-' e-(<VK B^)P^ . (28) 

and making use the definition of an average as 

— I dp M 

P^ = , 

Po 

we obtain the practical relation between mean pseudo-particle size and the IL density 

{n p i'^'f^ +C B'^)~p^ c^v , (29) 

where v = ^^-^ — . Apparently, it results in a well-known form of pseudo-particle size distribution 

M = C^J^^'^K'qcdp'-' ^ . (30) 

Now Eq.()29|) allows us to formulate the criterion we are interested in. It looks like C-B^ ^ n /? 
and for the IL parameters mentioned above it is 5 ^ 400 MeV. 
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Figure 3: The IL density as a function of applied external field. 



Dealing with Eq.(j2Zj) and Eq. ^Ujl the generating functional estimate Eq. ^Hjl may be presented 



as 



dx n < In 



n 



Y> I D[B] e-^(^) e-^ 



l-^+^^-ln[^(.)c^.^^^^ 



(31) 



.In^ 



Making use of the relation Eq. ()29p it is not difficult to find the maximum of functional Eq. ()31|) in 

dF 
dp 



the IL parameters at the fixed B value as a solution of transcendental equation = 0). As an 



information we give the simple expression of its derivative in n 

1 



n 



+ 



Fig.2 and Fig. 3 demonstrate the solutions for p and n at A'^c = 3 and Nf = 2 as the functions of 
field B. Fig. 4 shows the plot of free energy density f / ^qcd where F = J dx f and convinces IL is 
steady as to an impact of external field. At strong external field the IL parameters are given by the 
following asymptotic formulae 



P' 



n 



C B' 



1 - 



B^ 



2Nc 



T{u) Cn^ /?2Ar. N,buf3e 



2\u 



2\u 



e B' 
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This regime starts somewhere around BAqqj^ ~ 10 at all the plots given. 
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Figure 4: Free energy density as a function of external field B. 

Thus, the effective action for the B field is given by the following nonlinear functional 

'G%{B) G1^{B) 



S, 



eff 



dx 



4 



+ f[B] 



(32) 



This functional makes possible to calculate the external field as a function of x and IL parameters 
p[B] and n[B]. 

To get any estimate of the IL feedback on the presence of external field could be very practical 
for instanton liquid model. If so let us try to extract such an estimate from very simple example. 
Now we will search the minimum of effective action resolving the following boundary value problem 



ArB 
B 



df[B] 



dB 

r=ro Pi^) 



(33) 



Vj--B |r=ro 



e 1 

Air 



The source intensity here is controlled by e, and parameter tq sets a radius of colour ball which we 
take as ~ O.lp (albeit it is unessential) in order to avoid the difficulties in resolving the singular 
boundary value problem of Eq. (jH!?|l . The solution could be accomplished numerically probing such 
values of potential p{e) which provide with the solution going to zero magnitude at large values of r. 

The IL density as a function of r is plotted in Fig. 5 for ten various quantities of intensity. The 
extreme right hand side line corresponds to e/47r = 1 and the extreme left hand side corresponds to 
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r A 

Figure 5: The IL density a a function of r for different values of intensity. Extreme right hand side 
line corresponds to e/47r = 1. Going to the left corresponds to changing e/47r with a pace of 0.1 up 
to e/47r = 0.1 what corresponds to extreme left hand side hne. 

e/47r = 0.1. The same quantity of spacing corresponds to the hnes running to the right with intensity 
increasing. As it was expected the sohition has the Yukawa hkc behaviour which is weU seen in Fig. 6 
where ln(i?r) is plotted as a function of r for four various values of intensity with the pace of 0.1 and 
e/47r = 1 for the upper line. Fitting it with the linear function gives the estimate of screening radius 
which looks as follows 

Rd ~ (1.24 Aqcd)-' , 

Amazingly, this resuhs remains practicaUy unchanged for the whole interval of the intensities from 
e/47r = 0.1 to e/47r = 1 and implies that such a parameter characterizes (at least in this interval 
of values) the screening properties of IL itself. In a context of the model it looks like rather soft 
scale for the screening radius and might be taken as another confirmation of adiabatic approximation 
relevance for the Coulomb external field. Besides it hints the instanton vacuum could provide the 
significant energy loss of a parton at the later time of collision process contributing (and may be 
essentially) to jet quenching. 

Another point which is not studied here but worth of mentioning concerns the manifestation of 
instanton ensemble screening properties. It turns out rather unexpected because saturating instanton 
configuration is randomly oriented in colour space both with external field and with it switched 
off. Actually, anisotropy in colour space is routed (and is playing a role analogous to distribution 
function of colour charge) in the gluon field action of corresponding exponent, namely in the 'mixed 
component of gluon field. It is intuitively clear other instanton-like solutions (sensitive to the presence 
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of external colour field) could be even more adequate configurations. Such solutions were investigated 
in Ref. ini and was demonstrated that suitable scale and proper configurations affected by external 
field (crumpled instantons) could appear, indeed. 
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Figure 6: ln{Br) as a function of r for four various solutions. The upper line corresponds to 
e/47r = 1. Going down the lines correspond to decreasing e/An with spacing 0.1. 



Eventually let us comment on how it is essential that we are dealing with singular (anti-)instanton 
ensemble as a saturating configuration. Apparently, the screening properties of effective Lagrangian 
for external field B could be provided by any stochastic configuration of small characteristic size. 
The assumption of superposition ansatz validity occurs crucial to have all the leading contributions 
coming from the 'mixed' (repulsive) component of G{A, B) again. Another solution of the problem 
may appear, of course, in the quantum approach but this discussion is out of this paper scope. 
Studying the pseudo-particle behaviour while inside (anti-)instanton medium (n ^ 0) one could 
explore the interrelation of two mechanisms (the repulsive interaction and freezing the coupling 
constant out Ref. I7j) of fixing instanton size. 

The paper was supported by INTAS-04-84-398, NATO PDD(CP)-NUKR980668. The authors 
are also indebted to Barbara Jacak and Tetsuo Hatsuda for stimulating discussions. 
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